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Abstract. 

We investigate the semiclassical dynamics of massless Dirac fermions in 2+1 
dimensions in the presence of external electromagnetic fields. By generalizing the 
a matrices by two generators of the SU(2) group in the (25 + 1) dimensional 
representation and doing a certain scaling, we formulate a S — > oo limit where 
the orbital and the spinor degrees become classical. We solve for the classical 
trajectories for a free particle on a cylinder and a particle in a constant magnetic 
field. We compare the semiclassical spectrum, obtained by Bohr-Sommerfeld 
quantization with the exact quantum spectrum for low values of S. For the 
free particle, the semiclassical spectrum is exact. For the particle in a constant 
magnetic field, the semiclassical spectrum reproduces all the qualitative features 
of the exact quantum spectrum at all S. The quantitative fit for S = 1/2 is 
reasonably good. 



PACS numbers: 73.22.-f, 73.22.Pr, 03.65.Sq 

1. Introduction 

The recent developments in discovery of graphene p] and topological insulators [2] has 
aroused interest in the physics of Dirac fermions in two and three spatial dimensions. 
In device applications, it is often necessary to consider the system in the presence of 
external, spatially varying fields [3] 0]. While the Dirac equation generally can be 
solved numerically in these cases, the semiclassical limit of a quantum system often 
provides useful physical insights into the dynamics. Another interesting aspect of 
Dirac fermions is the strong coupling between orbital and spinor dynamics leading to 
the well known phenomenon of Zitterbewegung. This effect in solid-state systems, in 
particular graphene, has been a topic of recent interest [SJ [7] • 

The 2+1 dimensional Dirac equation can be looked upon as the analogue of the 
physical 3+1 dimensional Dirac equation. The Hilbert space of Dirac spinors are then 
projective unitary irreducible representations of the 2+1 Lorentz group [8]. In the 
condensed matter context mentioned above, there is no Lorentz invariance. The Dirac 
equation is valid only in a special frame, namely the rest frame of the material. The 
physical significance of the two components of the Dirac spinor depends on the details 
of the material as we describe below for graphene and topological insulators. 

The Dirac equation which describes the charge carriers in graphene is a 
consequence of graphene's honeycomb lattice crystal structure which consists of two 
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inequivalent triangular sublattices. The two bands arise due to the hopping between 
the two sublattices and touch each other at two points, the so called K and K points 
of the Brillouin zone. The energy dispersion is linear near these two points. The 
dynamics of the low energy quasi-particles and quasi-holes are effectively described 
by the massless 2+1 dimensional Dirac equation with the role of c played by the 
Fermi velocity vf- There are two species of Dirac fermions, one corresponding to the 
excitations near the K point and the other near the K point. For each of these, the 
two components of the Dirac spinor represent the probability amplitudes at the two 
sublattices. While the electrons in graphene are confined to two dimensions, they have 
a physical spin degree of freedom. This manifests in the effective Dirac theory as an 
internal degree of freedom. There are thus a total of four species of Dirac fermions in 
graphene. 

The 2+1 dimensional Dirac equation also describes the low energy quasi-particles 
at the surface of 3-dimensional topological insulators [H| ■ Examples of 3-dimensional 
topological insulators include Bi 2 Se3, Bi 2 Te3, Sb2Te3. While in graphene the Dirac 
equation arises from the honeycomb lattice structure, in topological insulators it comes 
from the strong spin-orbit coupling. Thus unlike graphene where the two components 
of the Dirac spinor represent the wave functions on the two sublattices, in topological 
insulators they represent the wave functions of the two spin projections. 

The earliest approach towards semiclassical analysis of a Dirac particle in an 
electromagnetic field was by W.Pauli[10 . He determined the phase of a WKB spinor 
for relativistic point particles for some special cases and argued that in the semi- 
classical limit, the translational motion is independent of the spin-degree of freedom. 
This conclusion was criticized by de BrogliejTT] with the remark that electromagnetic 
moments are expected to influence the trajectories in the semi-classical limit. The 
issue was latter addressed and clarified by Rubinow and Keller [12]. They pointed out 
that the moments of an electron are proportional to h so that in leading order as 
H the influence of the spin of the particle on the trajectories vanishes. Thus in 
a formalism when the spin is also taken to be large we can expect it to couple to the 
orbital dynamics. We analyze such a formalism in this paper. 

In their study of multi-component wavefunction Yabana and Horiuchi [T3] 
pointed out the necessity of including the geometric phase in the Bohr-Sommcrfcld 
quantization condition. However in their approach and many of the previous 
approaches to the semi-classical analysis of systems with many component wave 
functions, [HI [HI [16] the spin and orbital degrees of freedom are treated on a 
different footing, with the orbital degree of freedom being treated semi-classically and 
spin quantum mechanically. The orbital and spin degrees of freedom were combined 
together for the first time by Pletyukhov et. aZ[18] in a joint (extended) phase space 
where the description of a spin by continuous variables is achieved by using the basis 
of coherent states. The partition function for a system with spin-orbit coupling is then 
written in the path integral representation using coherent states. The semi-classical 
propagator or trace formula is evaluated for the system applying stationary phase 
approximation [17] . 

In this paper we follow the extended phase space approach and present a 
systematic formalism for the semi-classical analysis for the massless Dirac fermion 
in 2+1 dimensions in presence of external electromagnetic fields. We use the fact 
that, in 2+1 dimensions, the 3 Dirac matrices are elements of the two dimensional 
representation of the SU(2) algebra to generalise them to be the corresponding 
elements in the 2S +1 dimensional algebra. We then establish, using the path integral 
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formalism that there is a systematic stationary phase approximation with 1/5 being 
the small parameter. The classical equations of motion couples orbital and the spinor 
degrees of freedom. For cases of free particle on a cylinder and particle in a constant 
magnetic held closed orbit solutions are found and quantized using Bohr-Sommerfeld 
quantization. The classical solutions for the free particle are oscillatory rather than 
straight lines thus illustrating the phenomenon of zitterbewegung. Comparison with 
the exact quantum mechanical results shows that all the qualitative features of the 
spectrum are reproduced by the semi-classical analysis. The quantitative agreement 
is good even for small values of 5, in particular 5 = |. 

The organization of the rest of this paper is as follows. In section [2j we give the 
formulation of the path integral description of the propagator using coherent states and 
formulate a systematic semiclassical expansion in 1/5. We write down the classical 
equations of motion and the quantization condition for the closed classical orbits. In 
section [31 the Dirac fermion on a cylinder in absence of any external field is studied. 
The particle in a constant external magnetic field is studied in section 2] The semi- 
classical and exact spectrum are compared in section [5] The results are summarized 
in section [5] 

2. The Semi-Classical Expansion 

We begin by outlining the general formalism to analyze the spectrum semi-classically. 
The hamiltonian for a single massless 5 = i Dirac particle in 2 spatial dimensions, in 
the presence of external magnetic and spatially varying fields V(x) is given by: 

U = vpa.il + eV(x) (2.1) 

where, Ilj = pi — eAi is the mechanical momentum, vp is the Fermi velocity. The 
a matrices can be represented by the Pauli matrices and the wavefunctions by a two 
component spinor. 

The hamiltonian in eqn. |2.l| is generalized by replacing the a x and a v matrices 
by two of the generators of the SU (2) group in the (25 + 1) dimensional representation 
and written as: 

H = -^rS.II + eV(x) (2.2) 

We note that while the system is rotationally invariant at all 5, it is not Lorentz 
invariant for 5 ^ 1/2. We will explicitly see the consequences of this in the spectrum 
in later sections. However for applications to condensed matter systems, Lorentz 
invariance is not crucial. We resort to a path integral formalism to describe the 
time evolution of the system and find a systematic semiclassical parameter 5 for the 
combined orbital and spinor degree of freedom as described below. 

The path integral is constructed using coherent state representation. The coherent 
states which describe the system are direct product of the usual phase-space coherent 
states \x,p) and 5C/(2)-coherent states |n, /x) i.e 

\x,p;n,fi) = \x,p) <8> |n,/x) (2.3) 

where n.ft + /i 2 — 1. 

The path integral representation of time evolution of the system from time ti to time 
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t f is: 



(Xf,pf,n f ,fj,f |e R | Xi,pi]n h Hi) = / 2?[n(r, C)]2?[/x(r, C)] 



2?[x(t)]2?[p(i)]e— s" (2.4) 



f tf dx 

S[x(t);p(ty,n(T,0,Kr,{)}= J dt(f.— - U) - hSQ[n(r, C),/x(t, C)] (2.5) 
where 5 is the action with: 

x(^) = Xi] p{U) = p^, n(Ti) = n t ; ^(r,) = fa (2.6) 
= Pitf) = Pfi n( T f) = nf, K T f) = M/ (2-7) 
il[n(r, £), /z(r, £)] is the geometric phase encompassed during time evolution and is 
written as: 



f 



fi[n(r,C),M(T,C)] =J dr J^dC^d T rtd Q n\ i,j = 1,2 (2.8) 

fi(r; £) parameterizes the surface on the unit sphere in 2 dimensions defined by the 
curve n(r), and the geodesies from h{Ti) and h(rf) to the north pole (or some standard 
point on the unit sphere). For any arbitrary 5* it is convenient to write S = Sn. The 
system variables are scaled in terms of the characteristic length scale A and time scale 
lu of the system as: 

^Ts A - n; -A n - v '^ = ^ns v ^ (2 - 9) 

Using Eqn. |2.9j the action S can be written as: 
S[x'(t');p'(t');n(T 7 (),^{t, Q] 



= ns 



i i dx 

dt(p.— -m)-Q[n(T,0,Kr,0} 



(2.10) 



where, 



H = n.n +V (x) (2.11) 

Henceforth we will drop the subscript in the variables keeping in mind that they are 
scaled and hence are dimensionless. The propagator in ean. (12.41) can then be evaluated 
by a systematic stationary phase approximation with 1/5* as the expansion parameter. 
In this paper we will find closed orbit solutions to the classical equations of motion 
and quantize them using the Bohr-Sommerfeld quantization procedure. 

In the classical limit the system observables satisfy the following Poisson bracket 
relations: 

{x h p j } = 5 ij ; {ni,nj} = e^fi; {p,,m} = CjiUj {n h Xj}=0 
and the equations of motion are: 

Xi Til 

IT = e Jt jij - diV{xi), i = 1, 2 
hi = fieijUj, i,j = 1,2 

fi ^eijUiUj, i,j = 1,2 (2-12) 
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If the classical trajectories are closed then the action is quantized using Bohr- 
Sommerfeld quantization condition to find out the allowed orbits and subsequently 
quantize the spectrum. In terms of the canonical co-ordinate and the momentum, the 
quantization condition encloses an area in the phase space. For the Dirac particle in 
(2+1) dimensions, the quantization condition includes the contribution from geometric 
phase which is otherwise absent for quantization of Schrodinger particle, as follows: 



j> pidx + j> p 2 dy - Q[n, fi] = 2tt^7, 



m = 0, 1,2, • 



(2.13) 



where C denotes the closed orbits in the phase space. We use this method to analyze 
the spectrum for two solvable cases of (i) a free particle and (ii) a particle in a constant 
external magnetic field in following sections ([3]) and ((H) and compare with the exact 
solutions. 

3. Free Particle 



3.1. Semi- classical expansion 

Consider a massless Dirac particle on the surface of a cylinder of circumference L. In 
the absence of any spatially varying field V(x) and magnetic field, the hamiltonian for 
the system is: 

H /ree = n.p (3.1) 
We take A = L and uj = The action is: 



S[x(t);p(t);n,fj] = TiS 



dx 

dt\p x — - n^i] - Q[fi, fi] 



(3.2) 



The equations of motion are as follows: 

ft =o 

/".,/',• «,.?' = 1,2 
(jl = tijUiPj, i,j = 1,2 
We look into solutions for which p(t) = p®, p y (t) 



(3.3) 



and the solutions are: 



x(t) 

y(t) 
p(t) 

n x (t) 
n y (t) 



■ n x t + xq 

- —(u° cos uit + n° sin cut) 



n y cos ujt 
H° cos wt + n° sin uit 



H sin wt 

(3.4) 
where, uj = p x . 

The velocity and hence the classical trajectories thus oscillate in the direction 
transverse to the momentum illustrating the phenomenon of zitterbewegung. 
The geometrical phase traced out by the closed classical trajectories is: 

CI = 2tt(1 - cos(0 o )) (3.5) 
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COS (0 O ) =fi.P^ = n ° (3.6) 
From commensurate time-period for the linear motion and that for the orbital motion: 

2tt 1 



(3.7) 



The Bohr-Sommerfeld quantization condition following ean. (|2.13p for the free 
particle is: 

<j)p° x dx - 2tt(1 - cos 6>) = 2tt^ (3.8) 
where m is an integer. Substituting the value of cos 9 from (|3.6p . we have: 

p° x - 2tt(1 - n° x ) = 2vr^ (3.9) 

Quantizing each of the term individually, we get: 

p° x = 2nn/S (3.10) 

n° x =(mi/5 + l) (3.11) 

where n = 0, 1, 2, • • • , —25 < mi < 
The energy of the particle is given by: 

E = n x p% 

Tfl 

= 27rn( — ) where m = -S,...,S (3.12) 

Next we analyze the quantum mechanical spectrum of the free Dirac particle and 
compare with the spectrum obtained above. 



3.2. Quantum spectrum 

The hamiltonian in dimensionless units is given by: 

^ = ^S.p (3.13) 

This hamiltonian can be easily diagonalised. To compare with the semiclassical 
spectrum, we consider plane waves propagating in the x direction. The periodic 
boundary conditions imply that k x = 2im/L. S x has eigenvalues m = —S, . . . S. 

Tfi 

E = 2nn—, n = 0, 1,2- •• (3.14) 

Comparing equations (|3.12I) . (I3.14[) we find that the spectrum of the free particle 
obtained semi-classically is exactly identical to its quantum mechanical one. 

Also note that the quantum spectrum consists of 2S + 1 massless particles, 
travelling with speed tuvf- This illustrates our comment earlier that except at 
S = 1/2, the theory is not Lorentz invariant. 
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4. External magnetic field 

Semi- classical quantization 

In this section we study the massless Dirac particle in a constant magnetic field. The 
system is characterized by a length scale, the cyclotron radius l c and a time scale, the 
cyclotron frequency uj c which are: 

i - - vf (a-w 

lc — \ — w c = -— 4.1) 

V eB l c 

We choose A = l c and u> — lj c . The system variables are scaled in terms of the 
characteristic scales of the system following eqn. (|2.9l) . The system is described by: 



^mag 



n.n (4.2) 



where II related to the canonical momentum and the vector potential as IT = Pi + A4 . 
For a constant magnetic field, {II;, IL,} = e,j. We define another observable, the 
guiding center co-ordinate of the system Ri = Xi + CijUj, which satisfies: {Ri,Rj} = 
—£ij, {Ri, IX,} = & {Ri, rij} = 0. It then follows that Ri are constants of motion, 

{R i ,m mag } = ,i = l,2 (4.3) 

The canonical momentum pi and the canonical co-ordinate Xi can be written in terms 
of the guiding center co-ordinate and the kinetic momentum as follows: 

x = R x - U y ; y = R y + U x 

Px = {n x -R y )/2; Py = {n v +R x )/2 (4.4) 
The equations of motion of the system in terms of R , II ,n and n can be written as: 
R t =0 

hi = [ieijllj, i,j = 1,2 

a =e ij -n i II ! -, i,j = 1,2 (4.5) 

We choose Ri(t) = Ri. Since Ri commutes with the hamiltonian it represents the 
infinite degeneracy of each Landau level in the exact solution. The system has rotation 
symmetry in the x — y plane and for this closed system the angular momentum is 
conserved and is given as: 

^ = tijXiPj + n = ilTIT + n - ^RiRi (4.6) 

in terms of the guiding center co-ordinate and the kinetic momentum. 

The Bohr-Sommerfcld quantization condition in terms of II and Ri is given as: 

i j (eijUitlj - <, ,//,/>.,) alt - Q[n, /*] = 2tt^ (4.7) 

where Q,[n, fi] is given by eqn. (|3.5[ ). To solve for II, n and \x, we consider the following 
ansatze in which a) h = n x x + n y y parallel to II. b) n is perpendicular to the II. We 
look into each case separately in the following sections. 
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4-1.1. H || II For this case, we consider /j, = cos(0). The components of II are related 
to the azimuthal angle (f> in the x — y plane as: 

U x = ITcos</>; liy = rising (4.8) 

The equations of motion (|4.5j) in terms of 4> and becomes: 

sin00(f)=O; cos00(f) = O; tl(t) = 

sin 64> = -cos 011; 110(f) = 2 sin 6 (4.9) 
and the solutions are: 

0(f) = O (4.10) 

n(t) - n (4.11) 



0(f) = Vcos^o) « (4.12) 
where cos 9q = /io- The corresponding trajectories are: 
x(t) — R x = — IIo sin 0(f) 

y(i) - i? y = n cos 0(f) (4.13) 
where IIo is related to do as: 

sin^^ , /2<0o<7r (4.14) 

COS Uq 

The particle trajectories are circular orbits centered around the guiding center. 

For a given 9, the tip of the unit vector thus moves in closed orbits. The solid 
angle subtended by the closed path is given by eqn. (|3.5j) . We then use the Bohr- 
Sommerfeld quantization rule to pick the allowed orbits. The quantization condition 
is: 

II 2 , „ \ m , 

^+cos0 o j = — + 1, (4.15) 

where m = —2S, —2S + 1, • • • . With a change of variable n = m + 2S, where n is an 
integer and little algebra we get, 



cos O - I ((| - 1) - \J (| - I) 2 + 3j (4.16) 
from equations (|4.14[) and (|4.15[) . The energy of the particle is given by: 



£ = sin0oilo =± sin 2 O J 1 (4.17) 

y |cos0 o | 

1 n 2 r 2 

The angular momentum J z for this class of solution is 4f = + cos 0o — 

which using the Bohr-Sommerfcld quantization condition eqn. (j4.15|) . 

Jz (n \ R a 2 



S 

where Rq — RiRi. We demand (J z — S) to take integer values. This constraints 

Hi - k 
2 s 



= % where k is an integer. The quantum states corresponding to the different 
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values of k are degenerate, consistent with the infinite degeneracy of each Landau 
level. 

For large n limit, cos 8 ~ 2 n/s ' ener gy °f the particle in this limit is: 



e = n x ; ^ = n 9 -n x ^(sm(0) (4.22) 



E = ±y/2n/S (4.19) 

The energy spectrum for different values of S both for small values of n and in 
the asymptotic limit are plotted and discussed in section [5] 

4- 1.2. ft _L II For this ansatz, the components of II are related to the azimuthal 
angle <f> in the x — y plane as: 

U x = -rising; U y = II cos (4.20) 

so that the energy of the particle, given as : 

E = n x Ii y + n y Il x = (4.21) 

However, with this choice of co-ordinate axis, the variation of the polar angle 8 crosses 
singular points at poles. We choose a suitable co-ordinate frame to avoid the singular 
points. In the new frame the equations of motion are: 

tl x = — cos(#); ILy = sin(0) sin(</>) 
cos(0) ^ 
: sin(0) 

In the z — x plane (8 = ir/2) the above equations simplifies to: 

II = sin(0); = n (4.23) 
where II = II y . Combining the above equations we get, 

4> + sin ((f>) = (4.24) 

which maps to the pendulum equation. <f> maps to the angular co-ordinate or the 
amplitude, II to the momentum and cos cj> to the potential energy of the pendulum. 
The angular momentum 

^-cos0+^ (4.25) 

12 

maps to the total energy cos <fi + ^- of the pendulum which is a conserved quantity. 
The general solution to eqn. (|4.24p is Legendre's elliptic function of the first kind. II 
can also be subsequently computed from the values of <f> after subsequent substitution 
in eqn (|4~23l) . 

For J z < S the particle will execute libration. This involves a back-and- 
forth motion in real space and the phase space motion is contractable to a point 
in topological sense. The geometrical phase for this motion is zero. For J z > S the 
particle will undergo complete rotations in the z — x plane. The geometric phase 
picked up during one complete rotation is 2n. Since, = 0, the Bohr-Sommerfeld 
quantization condition for J z > S is: 
m 

n = 2tt— S = m (4.26) 
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where m is an integer. 

The zero energy modes for the regime J z < S occurs both for the integer and half 
integer S. However as the J z value increases and J z > S, the above Bohr-Sommerfeld 
quantization condition dictates that the zero modes, will appear only for the integer 
spins. 

To summarize, the semiclassical analysis predicts that there will be a zero energy 
state for every value of J z for integer S but for half odd-integer S, they will be zero 
energy states only for J z < S. 



4-2. Quantum spectrum 

In this section, we study the quantum mechanical spectrum of the system in a constant 
external magnetic field B. In the dimcnsionless units defined earlier, the hamiltonian 
for the system is: 

H = -=S ■ n (4.27) 

o 

with [H,IIj] = €ij. We define the operators a and as, 

a =-^(n s -in„) (4.28) 

a^ = ^=(U x +iU y ) (4.29) 

they satisfy the commutation relation [a,a^] = 1. With the standard notation, 
= S x ± iS y , the hamiltonian is: 

H = —= (S~a + + S+a) (4.30) 
<SV 2 

For S = 1/2 our model maps on to the strong coupling limit of Jaynes-Cummings 
model 12 . For general S, the model maps to the strong coupling limit of the Tavis- 
Cummings model (T9l [20] . 

The Tavis-Cummings model is exactly solved using the Bethe ansatz. However, we 
can deduce several qualitative features of the spectrum from the symmetry properties 
alone. The hamiltonian is rotationally invariant. The generator of rotations is, 

J z = a) a + S z (4.31) 

J z commutes with the hamiltonian and hence in the basis of the eigenstates of J 2 , 
the hamiltonian will be block diagonal. The eigenstates of J z can be constructed as 
follows. Let 

a)a\n, m) = \n, m) ,71 = 0,1,... (4.32) 
S z \n,m) —m\n,m) ) m = —S,...,S (4.33) 

we then have, 

J z \n-k,-S + k)=j\n-k,-S + k) ,j = n-S (4.34) 

For —S<j<S, there are n + 1 (= J + S + 1) values of k and thus the corresponding 
blocks of the hamiltonian will be J + S + 1 dimensional. For J > S the dimension of 
the blocks is 25 + 1. 
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For all values of S, there is an operator X = e ln " which anti-commutes with the 
hamiltonian. This implies that the eigenvalues of T-L have to come in pairs, namely 
if e is an eigenvalue then — e is also an eigenvalue. S commutes with J z and hence 
the eigenvalues of each block will have to come in pairs. Thus every block with odd 
dimension has to have at least one zero eigenvalue. 

The lowest value of j is j = —S. There is only one state in this block, |0, —S). It 
is easy to verify that this is always an zero energy eigenstate of H. Thus the blocks 
j = —S + 2n will all have atleast one zero eigenvalue for 2n < S. For j > S, the 
dimension of the blocks, as mentioned earlier, is 2S + 1. Thus for integer S, there 
is atleast one zero eigenvalue in each such block. For half-odd integer S, this is not 
necessary. 

The hamiltonian can also be diagonalised analytically in the j » S limit. In 
this limit, the hamiltonian reduces to 

H=^S y (l + o(l/n)) (4.35) 

where n = J + S. The energy eigenvalues of the 2S + 1 branches, in this limit are 
therefore, 

rrr 

E = V2n-£, m = -S,...,S (4.36) 

The energy eigenvalues can be numerically computed for low values of n. We 
have shown the spectrum for S = |, 1 and S = §, 5 in figure PQ. 

5. Discussion 

We now compare the results of the semiclassical quantization and the exact spectrum. 
In the free particle case the two spectra agree exactly. There is no exact match of the 
spectrum in a constant magnetic field and comparison is more interesting. 

We first consider the zero-energy states. The semi-classical analysis predicts 
that zero energy modes will occur for j < S. This is consistent with the results 
obtained from the exact analysis. The quantum spectrum has zero modes only the 
odd dimensional subspaces. This detail is not reproduced by the semiclassical analysis. 
However, we note that that the semiclassical regime in our formalism is the opposite 
limit, j > S. In this regime the semiclassical analysis predicts zero modes only for 
integer S and that is indeed bourne out by the exact spectrum. 

The exact spectrum has 25 + 1 branches for every j > S. However, we have 
only two branches of classical solutions. So which branch does it correspond to ? A 
comparison of the asymptotic results for the exact and the semiclassical asymptotic 
analysis, equations f|4.36[) and (|4.19|) shows that the classical solutions correspond to 
m = ±yo. Namely, should lie between ±m and ±(m + 1) where m < \S < rn + 1. 
Our numerical computations plotted in figure Q] shows that this is true even for small 
values of S. For S = 1/2, while the V2n dependence of the spectrum is reproduced, 
the semiclassical spectrum overestimates the energies by a factor « \[2. 

6. Results and conclusions 

To summarize our results, the semiclassical spectrum obtained by Bohr-Sommerfeld 
quantization reproduces the exact quantum spectrum for the case of a free particle on 
a cylinder. 
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Figure 1: (Color online) Plots showing the energy spectrum obtained from semi- 
classical analysis and the exact quantum spectrum. The curves with red square points 
are obtained from semi-classical analysis. The quantum spectrum is denoted by black 
triangles. Figs. ([Ta fl .([lb ]) . (|Ic ]) .([Id" l) are the plots for S =1/2, 1, 9/2 and 5 respectively. 
For the quantum spectrum, the energy branches corresponding to the different S z 
values for a given S are marked. To avoid clutter only the positive energy branches 
of the spectrum for high S values are plotted. In the large n limit, the semi-classical 
spectrum gives the y/S — th branch of the quantum spectrum. 

For the particle in a constant magnetic field, the semiclassical spectrum correctly 
predicts the occurrence of zero energy branches for integer S and their absence for 
half-odd integer S for j > S. The exact solutions have 25+1 branches in this regime. 
The classical solutions we have found yields a average of these branches. We note that 
the level spacing between these branches is ~ S a . Hence more detailed quantitative 
agreement may be obtained by doing the gaussian fluctuations about the classical 
solutions instead of Bohr-Sommerfeld quantization. 

We thank M.V.N. Murthy, Matthias Brack, Diptiman Sen and Bindu Bambah 
for useful discussions. M.M thanks M. Deka for discussions related to the numerical 
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